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Decay rates of unpolarized top quarks into longitudinally and transversally polarized W bosons 
are calculated to second order in the strong coupling constant a s . Including the finite bottom quark 
mass and electroweak effects, the Standard Model predictions for the W boson helicity fractions are 
T L = 0.687(5), F+ = 0.0017(1), and T- = 0.311(5). 
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There has been a continuing interest in the measure- 
ment of the helicity fractions of the W boson from top 
quark decays from the CDF collaboration pHD] and from 
the DO collaboration [cHlll] at the Tevatron at Fermilab. 

In the Standard Model (SM) the top quark decays pre- 
dominantly into a W + boson and a bottom quark. In- 
teresting observables, independent of the production rate 
that is difficult to predict precisely for a hadron collider, 
are the fractions of the three possible W helicities: Fl 
(longitudinal), J" + (transverse-plus) and T- (transverse- 
minus). In the leading order (LO) in the strong coupling 
constant a s (that is, without any gluon corrections), and 
in the limit of a massless bottom quark one has |T3 ] 



T L : T+ : T- = 



l + 2a; 2 







2x 2 



l + 2x 2 



(1) 



with Tj, + J-+ + T- = 1 and x = m\y/m t . Using mw = 
80.401(43) GeV [l| and m t = 172.8(1.3) GeV [3 we 



2 M 



Mr 



0.216(3) and T L : F+ : T- ~ 0.7 



get x 
0:0.3. 

The leading order decay t — > bW is a two-body pro- 
cess. With the V — A interaction, a massless b quark is 
left-handed, thus the W can only be left-handed or lon- 
gitudinal due to angular momentum conservation. One 
therefore has F + = provided no gluons are emitted. 

The above LO predictions are only marginally changed 
by the bottom mass. For a pole mass of m& =4.8 GeV 
one finds that the total rate T decreases by about a quar- 
ter per cent compared to the massless b limit. The helic- 
ity fraction Ft, slightly decreases while J-- increases, by 
about one per mil. The leakage into the transverse-plus 
fraction J r + is less than half per mil. Radiative correc- 
tions are a more important source of the transverse-plus 
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rate. However, as we shall see, when NLO and NNLO 
gluon radiation is included, F+ still does not exceed two 
per mil. Since only hard gluon emission can influence the 
helicity fractions, this smallness is a reliable prediction of 
the Standard Model. 

For this reason, the transverse-plus fraction J 7 + is a 
sensitive probe of New Physics effects such as a right 
chiral admixture to the SM current. The left and right 
chiral contributions do not interfere for m^ = leading 
to a quadratic dependence on the admixture parameter. 
The contribution of the right chiral contribution can be 
obtained from Eq. |T]) by exchanging J- + -H- J r _ whereas 
Tl remains unchanged. We mention that there are some 
indirect model dependent constraints on a possible right 
chiral admixture to the SM current from measurements 
of b ->• s + 7 decays p^Hl|. 

Let us summarize the theoretical prediction for the he- 
licity fractions. In addition to the O (a^) effects com- 
puted in this paper, we include the lower-order contri- 
bution (Born + O (a s )) [P9l - l2ll |. the leading electroweak 
corrections [22| and, for F + , the m& effect. The errors 
resulting from uncertainties in rrit,b,w and a s and an esti- 
mate of the higher-order effects are added in quadrature. 
We find 



T L = 0.687(5), 
T+ = 0.0017(1) 
T- = 0.311(5). 



(2) 



The relative errors for Tt, and J-- are small (of 0(1%)) 
and, for the largest part, result from the experimental 
error on the top mass. The error for F+ arises from 
uncertainty in a s and, to a lesser degree, in m^. Its 
absolute value is small but the relative error is large due 
to the fact that J-+ vanishes at LO for mj = 0. 

Various methods have been used by the CDF and DO 
collaborations to experimentally extract the helicity frac- 
tions from the top quark decay data (see the recent re- 
view [23|]). While previous analyses have performed two 



fits keeping one of the helicity fractions at its SM value 
more recent analyses measure the fractions .TV and J r + 
simultaneously in a two-dimensional fit p], [Tl|. Using 
such a model independent analysis the CDF collabora- 
tion quotes values of Fl = 0.88 ± 0.11(stat)± 0.06(syst) 
and F+ = -0.15 ± 0.07(stat) ± 0.06(syst) 7]. In a sim- 
ilar analysis the DO collaboration obtains Tl = 0.425 ± 
0.166(stat) ± 0.102(syst) and F+ = 0.119 ± 0.090(stat) ± 
0.053(syst) [11]. Both measurements are consistent with 
the SM predictions. 

The experimental errors on the helicity fraction mea- 
surements are still rather large but will be much reduced 
when larger data samples become available in the future 
from the Tevatron and from the LHC. Optimistically the 
measurement errors can eventually be reduced to below 
one per cent. For example, an early Monte Carlo (MC) 
study quotes measurement uncertainties of AJ 7 ^ = 0.007 
and A.^ = 0.003 for an integrated luminosity of 100 



fb _i at Tevatron II energies [24|. The corresponding 
event rates can easily be reached at the LHC within one 
year. A more recent MC study based on 10 fb _1 at the 
LHC quotes measurement uncertainties of AJ 7 ^ = 0.019, 
AJ"_ = 0.018 and AJ"+ = 0.0021 [l^. 

The improvements in the accuracy of the experimen- 
tal measurements have to be matched by corresponding 
advances in the theoretical sector. The NLO O (a s ) cor- 
rections to the helicity fractions were calculated in [19[ — 
l2l| . They lower Tl and increase T- by about one and 
two per cent, respectively, relative to their LO values. At 
NLO there is now a small contribution to the transverse- 
plus fraction J° + of 0.001. The corresponding NLO elec- 
troweak and finite width corrections were determined 
in [22j. They are smaller than the strong corrections 
and tend to cancel each other for both Fl and T- . 

It is desirable to improve the accuracy of the theo- 
retical predictions and to check the convergence of the 
perturbative series by computing the helicity fractions at 
NNLO. A first step in this direction was taken in [25|,[26[ x 
where the NNLO corrections to the total rate were found, 
exploiting the smallness of x = mw /mt- A series in pow- 
ers and logarithms of x was obtained and found to con- 
verge rapidly, so that its first few terms suffice. The aim 
of this paper is to use similar techniques to calculate the 
NNLO strong corrections to the three helicity fractions. 

We first determine the rate Tl of the top decay with 
longitudinally polarized W, replacing the full sum over 
W polarizations by a projector described below. The 
previous knowledge of the total rate is used to calculate 
the transverse rate Tt = L + + L_ from the difference 
Ft = r — L^. We use another projector to find the 
difference L + — L_. Finally, the helicity fractions Ti = 
Ti/T are determined. 

Our calculation follows the approach outlined in 
Ref. [201 . Using the optical theorem we compute the 




FIG. 1: Sample three-loop diagrams. Thick and thin lines 
denote top and bottom quarks, respectively. Wavy lines de- 
note W bosons and curly lines denote gluons. In the closed 
fermion loop all quark flavors have to be considered. 



decay width from the imaginary part of top-quark self- 
energy diagrams, 



in, 



■Im(S), 



(3) 



where £ denotes the one-particle irreducible self-energy 
diagrams. Sample diagrams are shown in Fig. Q] The 
unitary gauge is used for the W boson so that diagrams 
with Goldstone bosons are not needed. However, the 
i?£ gauge is used for the gluons with an arbitrary gauge 
parameter. The gauge-parameter dependence cancels in 
the final result. 

Since we set the mass of the bottom quark to zero, the 
integrals contain two scales, mt and mw- To reduce these 
integrals to single scale integrals, we use the method of 
expansion by regions (see, e.g., Ref. [28j|). In the present 
case, there are two regions to be considered. In the so- 
called hard region, the loop momenta are of the order 
of mt, while they are of order mw in the so-called soft 
region. The integrals become scaleless and vanish if a 
gluon momentum is soft. Thus, we are left with two 
contributions to each integral: one where all momenta 
are hard and one where only the IV-boson momentum 
is soft. For each contribution we construct appropriate 
expansions in the corresponding small quantities. The 
remaining single scale integrals are further reduced to so- 
called master integrals using Laporta's algorithm [29Tl30j |. 

Compared to the NNLO calculation of the total width, 
to get the partial rates with various W polarizations re- 
quires replacing the total rate projector 2 



T,p,V 



-<] 



[lis 



q^q u 



' w 



(4) 



by the longitudinal projector P^ or the transverse- 
plus/minus projectors P± / . The longitudinal projector 



1 An approximate value for the total rate was found in Ref. [27 



2 In the unitary gauge, the W boson propagator reads iP M "/(g 2 — 
m?„V 



reads |2C 



P£" 



(m 2 ^ -p-qq^){m 2 w p v ' -p-qq v ) 



"m 2 w m 2 \q\ 2 



(5) 



where p is the top quark momentum and q — (go; 0, 0, \q\) 
is the momentum of the W boson which propagates in the 
z-direction. The transverse projectors can be obtained 
with the help of the forward-backward projector [20J 



1 



mm 



ie^°Pp a q p . 



(6) 



with n € {1/2, 1}. d = A — 2e is the number of dimen- 
sions. Integrals of this type can be easily evaluated by 
performing the integrations over the time- like and space- 
like momentum components separately. 

For traces with an odd power of 75, we use the pre- 
scription of Ref. [32j and replace 



7^75 -> ^y£ M a,3,57 Q 7 / V • 



(10) 



The e tensor is stripped off and absorbed into the pro- 
jector. As a consequence the renormalization constant of 
the axial- vector current at the requisite order becomes 



One has P£" 



[W v -W^ ±Wp) /2. 

Partial helicity rates involve two technical challenges 
absent in the total rate calculation. First, there is an 
additional propagator-like structure l/|<f | n , n £ {1, 2} in 
Eqs. (|3J| and ©, and second, we have to deal with the 
presence of 75-odd traces in dimensional regularization. 
Our approach to both issues is outlined below. 

In the hard region, we express \q\ 2 through the prop- 



agator factor N = {p + qY 



mi 



2pq + q l as follows 
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[N 2 



(2p-g) 2 
2m 2 w N 
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m w 



Am^m^] 
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In Eq. (|7|) we use the fact that we are only interested in 
the imaginary part and that q 2 = m w on the cut. Now 
we can construct the desired expansions in myy/mt as 
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W\ 
1 



2m t y-r 

i=0 



2m 2 w N- 



Am 2 t m 2 w 



AN 2 



Am^ 
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i=0 



2m l w N 



Am t m w 



N 2 



(8) 



which we truncate at some order. Thus, the additional 
propagator-like structure from the projector is trans- 
formed into a scalar on-shell propagator with momentum 
p + q and mass m tl raised to arbitrary, integer powers. 
For the calculation of the polarized decays we need, next 
to the master integrals of Refs. [26|,|3l]]> twelve additional 
three-loop master integrals. 

In the soft region, we cannot perform an expansion of 
\q\, since |<f | 2 = q 2 — m^ and go is of order mw in the 
soft region. However, in this region the VF-boson loop 
factorizes. Therefore, we only have to replace the usual 
one-loop tadpole integrals with integrals of the type 



d d q 



(q 2 



<'W 
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Mr 



\n ' 



(9) 



Z A = l 



11 



1 



-f) l24 C ' C *-6 C ' T ' n 'h» (U) 



1 



where C F = (N 2 - 1)/(2JV C ) and C A = iV c are the 
Casimir operators of the fundamental and adjoint rep- 
resentation of SU(iV c ), respectively. For QCD we have 
N c = 3 and Tp = 1/2. rif denotes the number of quark 
flavors. Additionally, we have to include the finite renor- 
malization constant 

, 5 = 1 -> + (^) 2 (^_-^ 



+ y 6 C F T F n f 



(12) 



to restore the anti-commutativity of 75 . Both renormal- 
ization constants were determined at the three-loop level 
in Ref. [33|. 

A NLO check of the new methods used in this paper 
is afforded by comparing with the expanded form of the 
known NLO closed form results given in [19h21| . We 
found agreement up to O (x 16 ) . 

We present our results in terms of the reduced helicity 
rates T,- where 



Tt = 



G F m 3 t \V tb \ 
8V2V 



(13) 



with i 6 {L, +, — }. G F is Fermi's constant and Va 
is the element of the Cabibbo-Kobayashi-Maskawa ma- 
trix which governs transitions between bottom and top 
quarks. 

The analytical results of our calculation are too long 
to be presented here. Instead we present their numeri- 
cal values. In Table U we show successive terms in the 
power series expansion (in terms of [x n ] := (a;",x n lna;)) 
of the NNLO correction to the reduced helicity rates up 
to terms of order O ([a; 10 ]). Noteworthy is the fact that 
there are also odd powers in the expansion in [x] . These 
terms appear in the expansion of the parity-odd helicity 
structure function Tp 4 and thereby in T± (see Table HJ. 



p = (mt; 0,0,0) in the rest frame of the top quark such that 
p-q = rntgo- 



4 This follows the pattern in unpolarized and polarized top quark 
decays where the expansion of the five parity-even structure func- 
tions have n=even whereas the expansion of the five parity-odd 
structure functions have n=even/odd |2(J . 



The odd powers of x stem solely from the soft region 
of Tp. The leading x° contributions of 1^ and T are 
equal to each other. This is a consequence of the Gold- 
stone boson equivalence theorem. Between the [x 4 ] and 
[x e ] terms the power series expansion is somewhat erratic 
for r + , T_ and T. However, expanding up to O ([x 10 ]) 
Table Q] shows that one has sufficient numerical stability 
and precision for all three helicity rates and their sum. 
The contribution of the O ([a; 10 ]) term amounts to about 

0.01%, 1%, 0.06% and 0.03% of the total for f L ,f+,f_ 
and r, respectively. The convergence is slowest for T + . 
But then T + is numerically very small. 



TABLE I: Numerical values for O ([x n ] := (x n ,x n In a;)) terms 
in the ^-expansion of the NNLO corrections to the reduced 
partial and total helicity rates Ti 





Tl 


f+ 


f_ 
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\x°] 


- 1.958 -10 _a 








-1.958-10" 2 


[x 2 ] 


4.737- 10~ 3 


3.860 -10" 4 


-3.861 -10" 3 


1.262- 10 -3 


[x 4 ] 


6.710- 10" 4 


1.351 


10~ 4 


-9.917 


10~ 4 


-1.856- 10 -4 
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-5.339 
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5.339 


10" 4 
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-1.467- 10" 4 


1.186 


10~ 4 


4.878 


10~ 4 


4.597- 10 -4 


[x 7 ] 





7.696 


10~ 5 


-7.696 


10" 5 





[A 


- 1.702 -lO -8 


-2.333 


10~ 5 


-1.723 


10~ 5 


-5.758- 10" 5 


[A 





3.408 


10" 6 


-3.408 


10" 6 





[x 10 ] 


-1.274- 10 -6 


-1.844 


10~ 6 


-2.176 


10~ B 


-5.294- 10 -6 
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-1.434-10 -2 


1.610-10 -4 


-3.931 -10 -3 


-1.811 -10 -2 



In order to present our numerical results on the helic- 
ity fractions we define helicity fractions up to O (n) by 
writing (n = 0, 1,2 denote the contributions up to LO, 
NLO and NNLO, respectively) 



T, 



(n) 



U) 



T n r 



(14) 



where i = 



L,+,— . We further define the increments 



T. 



(n) 



T. 



(n-l) 



SF, 



(n) 



A^py 1 ' /^F> . We present our numerical results 



in the form J 7 , = T, 



(()) 



AT, 



j; (0) ^ 0, as J, 



-(0) 



pra+^r+sT 



(i) 



AJ, 



(2) 



r(2) 



and also, if 
For our nu- 



merical results we use a s (m t ) = 0.1073(24), which we 
obtained with the program RunDec [34] from the values 
a s {m z ) = 0.1176(20) and m z = 91.1876(21) GeV [1|. 
We find 

T L = 0.6978-0.0075-0.0023 
= 0.6978(1-0.0108-0.0033), 
+ 0.00103 + 0.00023, 
0.3022 + 0.0065 + 0.0021 

= 0.3022(1 + 0.0215 + 0.0070). (15) 

The results in Eq. (|15|) contain higher orders in a s from 
the expansion of the denominators in Eq. (|14[) . In order 
to maintain the constraint Tl + J~+ + F— = 1 , we prefer 
the unexpanded definition of helicity fractions (|14|l . 

The numbers in Eq. (1151) show the good convergence 
of the perturbative expansion, even though AT> /F> 



T- 



and the relative increments 



(for i = L, — ) is much smaller than AJ 7 ? ' / AF> ' . The 
NLO corrections to the helicity fractions are already close 
to the expected future experimental sensitivities and the 
NNLO corrections increase these by approximately a 
third. In particular, the NNLO calculation of the helic- 
ity fraction T+ remains at the order of 0.001. Should a 
measurement reveal a significantly larger value, it would 
be a clear signal of New Physics. 
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